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GPR55 is a newly de-orphanized cannabinoid receptor which belongs to the class A G- 
protein coupled receptors (GPCRs) family and binds constitutes of the plant, Cannabis sativa. It 
has been suggested that the manipulation of GPR55 may have a therapeutic potential in the 
treatment of inflammatory and neuropathic pain. 

The purpose of the present study was to refine the transmembrane helices (TMH) 
conformation in GPR55 that have significant sequence divergence from other class A GPCRs. 
The methods used were conformational memories (CM) to refine the transmembrane helices 
(TMH) of TMH2, TMH5, TMH6 and TMH7. The results of these calculations were used to 
modify the GPR55 computer model initially built in Reggio lab based on a rhodopsin template to 
generate refined inactive and active models of GPR55. 

The average proline kink angle and standard deviation for each set of conformational 
results generated by CM were measured using the Prokink program. A statistical analysis of the 
resultant face shift, wobble angle and bend angle of the helices containing proline was performed 
using the one sample t test and compared to the P-2 adrenergic receptor. 

The refined model of the inactive receptor of GPR55 obtained from the conformational memories 
is shown in figure 1 . In conclusion the results of conformational memories are consistent with the 
proposal of Ballesteros, that even though the overall structure of rhodopsin and other class A 
GPCRs may be very similar, there are localized regions where the structures of these receptors 
diverge. A significant range of conformational diversity could be generated by the presence of 
Pro-kinks and Cys/Ser/Thr residues. The results obtained should help to define the mechanism of 



drug receptor interaction relevant to cannabinoid physiological and pathophysiological functions 
including drug abuse. 
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CHAPTER I 
INTRODUCTION 



For decades marijuana remains one of the most widely used and abused drug 
around the world. It interacts with the brain, producing psychoactive and potentially 
therapeutic effects. Over the past decade it has becoming important to understand the 
mechanism whereby marijuana interacts in the human body. 

Marijuana is the buds and leaves of the Cannabis sativa plant. This plant contains 
more than 400 chemicals, including delta-9-tetrahydrocannabinol (THC), which is the 
main psychoactive chemical in this plant. THC is known to affect the short term memory, 
motor coordination, increases heart rate and levels of anxiety. 

It has being discovery the existence of naturally occurring marijuana-like 
substances in the human body. The discovery of specific receptors that are activated by 
smoking marijuana and the finding of endogenous cannabinoids, which can also activate 
this receptors have implicated the marijuana-cannabinoid research into mainstream 
science with important implications in human health and disease. 

Cannabinoid receptors are activated by a neurotransmitter called anandamide 
which belongs to a group of chemicals called cannabinoids, like THC. THC mimics the 
action of anandamide binding to the cannabinoid receptor and activates the neurons 
causing adverse effects on the mind and the body. 
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The cannabinoids have been suggested to play several roles in animals and 
humans, including the regulation of cell development and growth, modulation of immune 
function, nervous function, reproductive and feeding behavior, and in a number of 
biological functions from movement, memory, learning, appetite stimulation, and pain, to 
emotions, blood pressure, intraocular pressure, etc. The promiscuous action and 
distribution of cannabinoid receptors in most biological systems provides and extensive 
signaling capabilities for cross talk between different families of receptors, which may 
explain the numerous behavioral effects associated with smoking marijuana (1). 

At least two subtypes of cannabinoid receptor exist which have been cloned from 
animal or human source. These subtypes are termed CB1 and CB2. The CB1 receptor is 
present in the brain with highest density in the hippocampus, cerebellum and striatum. 
The CB2 receptor is found predominantly in the spleen and in the haemopoietic cells. 
The existence of this receptor provided the molecular basis for the immunosuppressive 
actions of marijuana. 

Cannabinoid receptors are members of the G-protein coupled receptor (GPCR) 
super family and mediate inhibition of adenylyl cyclase activity and reduce cyclic AMP 
levels. 

GPR55 

The goal of this research is to determine the structural features of the newly de- 
orphanized cannabinoid receptor GPR55 (2,3,4). This receptor belongs to the Class A G- 
protein-coupled receptor (GPCR) family and binds constituents of the plant, Cannabis 
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sativa. The CB1 and CB2 receptors were the first cannabinoid receptors discovered. 
There is biochemical, electrophysiological and behavioral data to suggest that even when 
the CB1 and CB2 receptors are knocked out (in CB1/CB2 double knockout mice), the 
endogenous cannabinoid anandamide still produces effects, suggesting the presence of 
additional cannabinoid receptor subtypes (5). 

The ability of GPR55 to recognize CB1 antagonist (SR141716A, Rimonabant) 
and act as agonist has been reported (2,6). However the aminoalkylindole WIN55, 212-2, 
which is a CB1/CB2 agonist, did not produce an agonist response at GPR55 indicating 
that this is a distinct receptor from CB 1(4,7,8). Palmitoylethanolamide (PEA) has been 
identified as a potent anti-inflammatory, anti-excitotoxic and anti-hyperalgesic 
compound, which can act as a endocannabinoid for the GPR55 (4,9,10). Thus GPR55 
recognizes many cannabinoids, but has a unique response profile differing from CB 1 and 
CB2. GPR55 exhibits low amino acid identity to CB1 (13.5%) or CB2 (14.4%) receptors. 
The most closely related proteins to GPR55 are GPR35 (27%), P2Y (29%), GPR23 
(30%), and CCR4 (23%) (11). 

Ryberg and co-workers have suggested that the role that GPR55 plays in health 
and disease and any potential therapeutic benefits of activation or blocking of this 
receptor are yet to be ascertained (12). More recently, the fact that that GPR55 signaling 
can influence the regulation of certain cytokines has led to the hypothesis that GPR55 
may contribute to the lack of inflammatory mechanical hyperalgesia in the GPR55 ~'~ 
(knockout) mice. Consequently, it has been suggested that the manipulation of GPR55 
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may have therapeutic potential in the treatment of inflammatory and neuropathic pain 
(13). 

Two receptors for anandamide were found in human endothelial cells, characterized as 
cannabinoid 1 receptor (CNR1) and G-protein-coupled receptor 55 (GPR55). Both 
receptors trigger distinct signaling pathways upon anandamide stimulation where the 
integrins are activated and the signaling cascade becomes promoted. It has been reported 
that the effect of anandamide on endothelial cells depends on the status of the integrin 
clustering (14). 

GPR55 is highly expressed in larger dorsal root ganglion neurons and upon 
activation by various cannabinoids (A 9 THC, the anandamide analog R-methanandamide 
and JWH015) increases intracellular calcium in these neurons. The signaling pathways 
involved in HEK293 cells expressing GPR55 which lead to an increase of calcium 
involves G q , G12, RhoA, actin, phospholipase C, and calcium released from IP3-R gates 
stores. These pathways are quite different from those used by CB1 and CB2 (15). 

Lysophosphatidylinositol (LPA) has been reported to be a possible endogenous 
ligand for GPR55 (see Figure 1). In HEK293 cells expressing GPR55, LPA produces a 
rapid Ca 2+ transient signal and stimulation of GTPyS(16). Pertwee and co-workers 
reported that both the CB1 receptor antagonist, AM251 and the CB1/CB2 agonist, 
A 9 THC are able to activate GPR55. In addition, the CB1 endogenous antagonist 
virodhamine and non-CBl/CB2 ethanolamides, palmitoylethanolamide and 
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oleoylethanolamide activate GPR55 by coupling with Gan (16). Information about 
additional GPR55 ligands is summarized in Table 1. 




Figure 1 . L-a-Lysophosphatidylinositol (Liver, Bovine-Sodium Salt) 

Pertwee and co-workers also compared the distribution of GPR55 and CB 1 
mRNA in mouse (16) (see figure 2). The results of the expression profile of GPR55 in a 
panel of mouse tissue using PCR showed the presence of GPR55 mRNA with high levels 
in the adrenals, parts of the gastrointestinal tract and the CNS. Similar to CB1, GPR55 
has a broad distribution in brain tissue but the levels are significantly lower than those for 
CB1(16). 
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Figure 2. mRNA expression levels of GPR55 and CB1 receptor in mouse tissues 
measured by quantitative PCR. Data are mean values + s.e.m. using tissue from eight 
(GPR55) or four mice (CB1). 
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Table 1. Comparative table of the different GPR55 ligands. 



The orphan receptor GPR55 is a novel 
cannabinoid receptor (E. Ryber et al, 
September 17 th , 2007) 


Identification of GPR55 as a 
lysophosphatdylinositol receptor (Saori 
Oka et al., August 24 th , 2007) 


The GPR55 Iigand L-a- 
lysophosphatidylinositol promotes RhoA- 
dependent Ca 2+ signaling and NAFT 
activation (Henstridge C., et al August 29 th , 
2008) 


Expression profiling 


Mouse tissues using PCR 
-adrenals 

-Part of intestinal tract 
-CNS 

But levels are lower than CB 1 




Radioligand binding assay 


-specific binding for CP55940 
-small specific binding for SR141716 
-no binding for WIN55,212-2 








2AG, anandamide, PEA, N- 
oleoylethanolamide, virhodamide, 
CP55940, HU210, WIN55,212-2, -A 9 
THC, abnormal cannabidiol, and 
SR141716 did not exert any effect 
-LPI induced a rapid phosphorilation of 
ERK cell expressing GPR55 








-LPI was capable of inducing receptor 
internalization from the plasma membrane to 
the intracellular vesicles 




-CP55940 increased GTPyS binding with 
EC 50 18nM 

-2-AG, noladin ether, PEA, virodhamine 


-LPI dose dependently stimulated the 
binding of 

t 35 S] GTPyS to the GPR55 expressing cell 





and OEA stimulated GTPyS binding 
-A 9 THC, HU210 AM251 stimulated 
binding 

-cannabidiol antagonize the agonist effect of 
CP55940 

-WIN55,212-2 and AM281 no functional 
activity 


membranes 




G protein coupling 


-No effect Pertussis toxin no Gi 
-No calcium signaling no Gq 
-G13 present GTPyS binding in a dose 
dependent manner 




-LPI mediated Ca 2 + signaling involving Gal 3 




30 nM of LPI provoked a Ca 2+ transient in 
GPR55-expressing cells 


-Exposure to LPI resulted in a concentration 
dependent induction of oscillating Ca 2+ -LPI 
evokes Ca 2+ released from intracellular stores 




i 




Treatment of cells with an increased amount of 
LPI resulted in a dose dependent induction of 
NFAT activity 




-Anandamide and 01602 but not 
WIN55,212-2 induced activation of rhoA, 
cdc42 and racl 

-Effect was blocked by cannabidiol 




-RhoA is a downstream effector of 

Gal2/Gal3 family 

-luM of LPI elicited RhoA activation 


Do cannabinoids promote GPR55- mediated Ca2+ signaling? (profile of activity of selected cannabinoids ligands in GPR55-HEK293 






-AEA and 2- AG were able to elevate 
intracellular Ca 2+ levels in a concentration 
range of 3-30 uM 

-but not difference in maximal response, 
duration and potency were observed compared 
to the control 





-CP55940 had no effect on intracellular Ca 2+ 
levels 

-AM25 1 was able to induce sustained Ca 2+ 
response but it was less potent than LPI 




-vascular tone -controlling receptor 
-inflammatory activities 


-LPI is an acidic lysophospholipid 
-LPI is a possible degradation product of 
phosphatidylinositol 

-LPI reduced the response induced by the 
subsequent stimulation with LPI, 
suggesting desensitization in LPI 
pretreated cells 

-levels of LPI in patient with ovarian 
cancer was markedly higher 
-LPI may stimulate the growth and 
metastasis of ovarian cancer 


-Characteristic of GPR55 are unique 

a) no other GPR55 links the signaling 
pathways of Gal 3 and small G-TPases like 
RhoA 

b) with oscillatorv and prolonged Ca 2+ 
signaling 

c) and downstream NFAT activation which can 
regulate DNA transcription and gene 
expression 

-LPI can activate GPR55 

-but endocannabinoids AEA and 2- AG have 

not clear effect at equivalent concentrations 

-GPR55 knockout animals data suggest a 

potential role within the brain, neurons and 

immune cells regulating 

inflammatory/neuropathic pain 

-LPI is known to have mitogenic activity and 

may influence the growth of certain tumors 



GPCR Structure and Function 

GPCRs are transmembrane proteins that serve as a link to cellular signal 
transduction mechanisms. Due to the technical difficulties associated with crystallizing 
transmembrane proteins such as GPCRs, only a few GPCRs have been crystallized to 
date. These include rhodopsin (17,18,19), the beta-2-adrenergic receptor ([32- 
AR)(Cherezov et al., 2007; Rasmussen et al., 2007; Rosenbaum et al., 2007), pl-AR 
(Warne et al., 2008), and adenosine A2A receptor (Jaakola et al., 2008). Rhodopsin is the 
prototypical Class A GPCR and its crystal structure represents the dark-adapted (inactive 
state) (17,18,19). It has been used as the basis for homology modeling of other Class A 
GPCRs. The general topology for the GPCRs include: an extracellular N terminus seven 
transmembrane (TM) alpha helices arranged to form a closed bundle; loops that extend 
intra- and extracellularly connecting TM helices an intracellular C terminus that begins 
with a short helical segment (helix 8) oriented parallel to the membrane surface. 
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Figure 3. General topology of GPCR's. 

Ligand binding is thought to occur within the binding site crevice formed by the 
transmembrane helix (TMH) bundle, to extracellular loops, or to a combination of 
extracellular loop and binding site crevice residues. The homology models developed in 
the Reggio lab for the CB1 and CB2 receptors have been modified based on sequence 
dictated divergences between rhodopsin and each cannabinoid receptor subtype. These 
models have been very successful in identifying key interaction sites for classical 
cannabinoids (20), non-classical cannabinoids (21), aminoalkylindoles (22) and 
biarylpyrazoles (23). 

The primary structural information about GPCR activation comes from 
biophysical studies of rhodopsin and the p2-adrenergic receptor (P2-AR). These studies 
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indicate that the rotation of both TMH 3 and TMH 6, as well as a conformational change 
in TMH 6 (straightening in the CWXP hinge region of the TMH6) occurs upon GPCR 
activation (24,25,26,27,28,29). The time scale for the activation of the P2-AR by its 
diffusible ligand has been reported as seconds (25). This has been shown to be true in live 
cells, as well (29). This time scale is quite long compared to the typical length of 
molecular dynamics (MD) simulations (10- 100 ns). For this reason different models of 
the active states of CB1/CB2 (22,30) and GPR35/55 had been created based on 
documented changes that occur during the activation process. Light activation of 
rhodopsin does not require the breaking and forming of thousand of specific contacts 
within nanoseconds, only specific contacts restricting the inactive state are broken during 
activation. These changes are then transmitted through the entire membrane protein 
because of its dynamic plasticity. It has been reported that a key tryptophan residue in the 
ligand binding pocket, W6. 48(265) is restricted in rhodopsin. In the dark (inactive) state 
of rhodopsin, the beta-ionone ring of its covalently bound ligand, 1 1-cis-retinal is close 
to W6.48(265) of the CWXP motif on TMH F and helps constrain it in a Xi= g+ 
conformation (17,18,19). In the presence of light, this receptor activates and the beta- 
ionone ring moves away from TMH F and toward TMH D where it resides close to 
A4. 58(169) (31). During this process the constraint on W6.48(265) is released, making it 
possible for W6.48(265) to undergo a conformational change. This suggests that the 
conformation of W6.48 (265) when rhodopsin is inactive (R; Xi= g+) changes during 
activation to a trans conformation (%i g+ — > trans) (32). Another interaction that has been 
reported as important is W6.48/F3.36 which may act as the "toggle switch" for 
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CB1 activation, with W6.48 %i g+ /F3.36 Xi trans representing the inactive state (R) and 
W6.48 %i trans /F3.36 %i g+ representing the active state (R*) state of CB1 (33). It has 
been suggested that F3.36 may restrain W6.48 from moving to an activated state 
conformation (34). 

Table 2: Summary of GPCR Activation 



GPCR ACTIVATION 


Inactive (R) State 


Active (R*) State 


W6.48xling+ (-60°) 
N-Ca-Cp-X 


W6.48 xl in trans (180°) 
N-Ca-Cp-X 


"Ionic lock": 
Rho: Arg 3.50 & Glu 6.30 


"Ionic lock" broken: 
TMH3/6 Rotate 


Toggle Switch (p 2-AR) 
C6.47/W6.48/F6.52 
trans g+ g+ 


Toggle Switch (p 2-AR) 
C6.47/W6.48/F6.52 
g+ trans trans 



To build the activated state model of GPR55 we used the biophysical literature on 
the R to R* transition in Class A GPCRs, such as rhodopsin and p2-adrenergic receptor. 
This literature suggests that activation is accompanied by a straightening of the TMH 6 
and rotation of TMH 3 and TMH 6. The identity of the residue 6.49 can have a profound 
effect on the flexibility of TMH 6 (35) GPR55 like rhodopsin has a Leu at position 6.49. 
To generate an activated model of GPR55 we selected a conformer for which the proline 
kink angle is moderated as compared to the TMH 6 conformation in the rhodopsin crystal 
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structure (19). TMH 6 was also rotated so that Cys 6.47 becomes accessible from inside 
the binding site crevice (27). TMH 3 was rotated so that residue 3.41 changes 
environments (26). We expect that the most significant changes in the R to R* transition 
will be seen in the binding pocket of GPR55 due to the change of residue accessibilities 
on TMH6. The groove residues V6.43/V6.46 will rotate to the TMH6/7 interface, S6.47 
will become accessible to the ligand binding pocket. And H6.52 and F6.59 will rotate out 
of the binding pocket. 

In the present research, the amino acid numbering sequence to be used is the one 
proposed by Ballesteros and Weinstein (36) In this numbering system, the most highly 
conserved residue in each transmembrane helix (TMH) is assigned a locant of 0.50. This 
number is preceded by the TMH number and may be followed in parentheses by the 
sequence number. All other residues in a TMH are numbered relative to this residue. The 
sequence numbers used here are from human GPR55 (37). In the Ballesteros-Weinstein 
numbering system, for example, the most highly conserved residue in TMH 2 of the 
GPR55 receptor is D2. 50(70). The residue that immediately precedes it is called 
F2. 49(69). Residues in the intracellular extension of TMH 7 are numbered as if they are 
part of TMH7 following the literature precedent set by Prioleau and co-workers (38). 
Among the highly conserved residues typically used in sequence alignments with 
rhodopsin, GPR55 has the conserved patterns in TMH1, 2, 4 and 5 (i.e. N1.50, D2.50, 
W4.50, and P5.50). The fact that CB1 and CB2 lack the highly conserved P5.50 while 
GPR55 has this motif, suggests a structural divergence in TMH5 between CB1/CB2 and 
GPR55. GPR55 has the substitution DRF for the conserved TMH3 E/DRY motif in 
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CB1/CB2. In TMH6, the highly conserved CWXP motif found in rhodopsin and 
CB1/CB2 is substituted with SFLP in GPR55. The greatest divergence from the 
rhodopsin sequence and from CB1 and CB2 appears in THM7 because the highly 
conserved NPXXY motif is replaced with DVFCY (GPR55). Like rhodopsin, GPR55 has 
an F in the intracellular extension of TMH7 (called Hx8) at position 7.60. In rhodopsin, 
there is an aromatic interaction between Y7.53 (306) and F7.60 (313), which has been 
proposed to provide structural constraints that rearrange in response to 
photoisomerization (39). In GPR55, the analogous relationship between Y7.53 and F7.60 
can be established, but this interaction is not possible in CB1/CB2 because at position 
7.60 there is Leu in CB1 and He in CB2. 

GPR55 also potentially has another significant similarity to rhodopsin in its EC-2 
loop structure. In rhodopsin, the EC-2 loop dips down into the binding pocket to form a 
disulfide bridge between an EC-2 Cys residue and C3.25. It then loops back over itself to 
make its connection with the top of TMH5. GPR55 has a Cys at 3.25 and a Cys residue in 
the EC-2 loop that could potentially form a disulfide bond. It is likely that EC-2 loop 
structure of GPR55 will differ from that of CB1 and CB2 as these latter receptors lack the 
Cys at position 3.25. 

The initial model for GPR55 built in the Reggio lab was based upon a rhodopsin 
template. Work with CB1 and CB2 has shown that although the rhodopsin crystal 
structure is a good starting point for cannabinoid receptor modeling, there are sequence 
dictated differences between rhodopsin and each of these receptors that can result in local 
divergence from the rhodopsin structure. These can be manifested by the absence of a 
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key motif, by divergent locations of proline or glycine residues, as well as by the 
presence of serines or threonines in key positions that can exert subtle, but important 
effects upon helix conformation. 
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CHAPTER II 
OBJECTIVE 

The objective of my thesis is to use the biased Monte Carlo/simulated annealing 
method, Conformational Memories (CM) to refine the transmembrane helix (TMH) 
conformations of TMHs in GPR55 that have significant sequence divergences from other 
Class A GPCRs. These include TMH2, TMH5, TMH6 and TMH7 of GPR55. The result 
of these calculations will be used to modify the GPR55 computer model initially built in 
the Reggio lab based on a rhodopsin template to generate refined R and R* models for 
GPR55. 
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CHAPTER III 
METHODS 



Molecular Modeling 

Definition of torsion angles 




Figure 4. Torsion angels 

Dihedral angle cp : rotation about the bond between N and Ca bond, involves 
C(0)-N- Ca-C(O) 

Dihedral angle y : rotation about the bond between Ca -C(O) bond involves N- 
Ca-C(0)-N 

Dihedral angle w : rotation about the peptide bond C(0)-N bond involves Cal- 
C(0)-N- Ca2 
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(0 angle tends to be planar (0° - cis, or 180 ° - trans) due to derealization of 
carbonyl pi electrons and nitrogen lone pair 

cp and \\i are flexible, therefore rotation occurs here 

Side chain conformation 

The Greek alphabet is used to named the side chain atoms of the amino acids 
according to this scheme. 




C 





epsilon 



C delta 




C 



gamma 




c 



alpha 



lysine 



Figure 5. Side chain conformations. 
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The side chain torsion angles are named cl(chil), c2(chi2), c3 (chi3), etc., 



gauch& conformation 




"CO" group 



Figure 6. Gauche plus conformation 



The most abundant conformation is gauche(+) in which the gamma side chain atom is 
opposite to the residue's main chain carbonyl group when viewed along the CP-Ca bond. 




trans conformation 



Figure 7. Trans conformation. 
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The second most abundant conformation is trans in which the side chain gamma atom is 
opposite the main chain nitrogen. 



gauche ~ conformation 



Chi1 

= 60 
deqrees 




'CO' group 



Figure 8. Gauche minus conformation 



The least abundant conformation is gauche(-). Occurs when the side chain is 
opposite the hydrogen substituent on the C a atom. This conformation is unstable 
because the gamma atom is in close contact with the main chain CO and NH groups. The 
gauche(-) conformation is occasionally adopted by serine or threonine residues in a helix 
where the steric hindrance is offset by a hydrogen bond between the gamma oxygen atom 
and the main chain. 

Models of TMH 2, TMH 5, TMH 6 and TMH 7 for GPR55 

Low free energy helix conformations for GPR55 TMH2, TMH5, TMH6 and 
TMH7 were created using the conformational memories (CM) technique to explore the 
sequence dictated conformations of these key helices. In most cases, the helix was built 
using the standard § (-63°) and \\f (-41.6°) backbone dihedrals for TM helices (Ballesteros 
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and Weinstein 1995) (36). Backbone phi and psi torsion angles were allowed to vary by + 
10°, with a larger variation (+ 50°) used in more flexible regions of the helix, i.e. those 
containing known helix bending residues, such as proline, glycine, serine or threonine 
(40). Side chain torsions were allowed to vary + 180°, unless otherwise specified. The 
protocol for bond angle variation will be discussed below. 

TMH2 Flexible Region: TMH2 in Rho does not have a proline, however it does 
have helix distortion from G2.56 to T2.59 (GGFT). GPR55 lacks this GGFT sequence, 
but it has a proline at 2.58. It is possible that this proline is a structural mimic of GGFT. 
In CM calculations, the region from V2.54 to P2.58 was considered to be the flexible 
region in GPR55 TMH2. S2.56 near P2.58 was set to be in g minus xl exploring this 
region and this helix will be build based on the backbone of the P2 adrenergic receptor 
conformation for the proline kink. 

TMH5: In Rhodopsin and the P2 adrenergic receptor, there is a single proline 
residue in TMH5 at position 5.50, which causes a deviation from normal alpha- helicity in 
the region, 5.46-5.50. TMH5 in GPR55 has two proline residues (P5.41 and P5.50) which 
could cause distortions from normal alpha-helicity. GPR55 TMH5 was built based on the 
backbone of the P2 adrenergic receptor in order to capture the helix distortion around 
P5.50. The region from K5.37 to P5.41 was then considered flexible in CM calculations. 

TMH6: TMH6 in GPR55 does not have the Class A GPCR highly conserved 
CWXP motif. However, it has the conservative substitutions of this motif: SFLP. Since 
TMH6 has been implicated in conformational changes during receptor activation 
(24,25,26,27,28,29), it was very important to study the allowed conformations of GPR55 
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TMH6. Results from these studies were used to build not only the inactive state model of 
GPR55, but also the activated state model. CM calculations for TMH6 considered the 
region from V6.46 to P6.50 to be flexible in GPR55. The residue serine 6.47 was run 
under two conditions. In the first condition, the S6.47 %l was not restricted. In the 
second condition, the S6.47 %l was held in trans position due to its position and possible 
structural role found by Juan Ballesteros (36). 

TMH7: TMH7 in Class A GPCRs usually has the conserved sequence motif, 
NPXXY, which not only influences the conformation of TMH7, but also places Y7.53 in 
the correct position to interact with F7.60 on the intracellular extension of TMH7, Hx8. 
GPR55 lacks the NPXXY motif in TMH7. Instead, this region of the sequence is DVFCY 
in GPR55. The CM protocol for CB1 in TMH7 calls for the region from S7.46 to P7.50 
to be considered the flexible region of TMH7. To parallel this, the span from C7.46 to 
V7.50 (GPR55) was considered to be the flexible region. 

Conformational Memories (CM) 

In order to explore the different conformation of the helices TM2, TM5, TM6 and 

TM7 the conformational memories (CM) method was used. This method employs 

multiple Monte Carlo/simulated annealing random walks and the CHARMM 32 force 

field (41). The CM method is widely used because it converges in a very practical 

number of steps and it is capable of overcoming energy barriers efficiently. The CM 

method can fully characterize the conformational properties of a helix by the free energy 

of each of the conformations that the helix can adopt, which includes not only the 

intrinsic energy of each conformational state, but also the probability that the helix will 
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adopt each particular conformation relative to all other ones accessible in an equilibrated 
thermodynamic environment. The calculations are performed in two phases. In the first 
phase repeated runs of Monte Carlo/simulated annealing are carried out to map the entire 
conformational space of the helix. The dihedral and bond angle space in which 
conformations were accepted are used to create "memories". In the second phase, new 
Monte Carlo/simulated annealing runs are performed only in the populated regions 
("memories") identified in the first phase of the calculation. All calculations are run using 
a distance dependent dielectric. 
The two phases are described in more detail below: 

Exploratory phase 

In the exploratory phase of the CM method, a random walk is used to identify the 
region of conformational space most probable for each torsion angle and bond angle. The 
initial temperature for each run is 3000 K with 50,000 Monte Carlo steps applied to each 
torsion or bond angle variation with cooling in 18 steps to a final temperature of 310 K. 
Each step consists of varying two dihedrals angles and one bond angle chosen at random 
from the entire set of variable angles. The torsion angles and bond angles are randomly 
picked at each temperature and each move is accepted or rejected using the Metropolis 
criterion (30). Accepted conformations in the Exploratory Phase are used to create 
"memories" of torsion angles and bond angles that were accepted. This information 
provides a map of the accessible conformational space of each TMH as a function of 
temperature. 
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S18 N-CA-CB-OG ChM 




5.000 scale: 1.00000. 1.00000 

Figure 9. Map of the accessible conformational space as a function of temperature. CM 
Map of S4.54 Chil Pre-biased (Exploratory) Phase. 

Biased annealing phase 

In the second phase of the CM calculations, the only torsion angles and bond 
angles moves attempted are those that would keep the angle in the "populated 
conformational space" mapped in the exploratory phase. The biased annealing phase 
began at 749.4 K cooling to 310 K in 7 steps. Finally, 105 structures were output at 310 
K for each TMH. 

Metropolis criterion 

In the metropolis criterion or algorithm the first conformation is randomly 
generated and each point in the construction of the conformation a move is attempted to 
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the current conformation. If the local chain conformation is not compatible with the 
attempted move then the movement is rejected immediately. If the conditions are 
satisfied then the metropolis criterion is applied. The acceptance probability (P) of a trial 
move is chosen to be 

P = min(l;exp (-0 «AE) 

where (3= 1/T. 

Variable bond angles 

The memories for bond angles are defined in the range 90° to 144° with 1.5° bin 
sizes, using pre-biased memories for all bond angles to avoid sampling high energy 
regions of bond angle space. The bonds angles present in aromatic rings and side chain 
bond angles that include non polar hydrogens are not varied. The other angles are 
allowed to vary in a range of + 8° around the initial value. Bond angles involving polar 
hydrogens or the C-S-C angles are allowed to vary up to + 15° around the initial value. 
The full set of bond angles and dihedrals angles of a molecule determines the molecular 
structure and using the internal coordinates of CHARMM, the minimum set of bond 
angles necessary to construct the coordinates of the system are determined. 

Evaluation of CM Output 

Helix Geometry. The average proline kink angle and standard deviation for each 
set of conformational results generated by CM were measured using the ProKink 
program (40). This program, which is part of Simulaid Conversion program was used to 
calculate the face shift, wobble angle, and bend angle of each helix. A statistical analysis 
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of the resultant face shift, wobble angle and bend angle for the 105 helices was performed 
using the statistical analysis package (GraphPad Software for one sample t test). 

Helix Choice for TMH Bundle. 

The 105 helices output for each run were superimposed on the appropriate helix 
in the beta-2 adrenergic receptor template. For TMHs 2,5, 7 and the TMH6 inactive state, 
output helices were superimposed from the intracellular side up to (but not including) the 
flexible region. For the TMH6 activated state, CM output was superimposed from the 
extracellular end up to (but not including) the flexible region. For all helices except 
TMH6, helices were chosen based upon their fit within the TMH bundle. The helix 
chosen for the TMH6 inactive state was chosen for fit and also for its ability to form an 
R3.50/Q6.30 interaction. The helix chosen for the activated state of TMH6 was chosen 
such that this interaction was broken. Additional criteria used for helix selection will be 
discussed in the Results section. 
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Figure 10. Presence of serine in a Helices 

Ballesteros et. al., documented the effects of Serine residues found in a g- or +60° 
where the hydrogen bonding to the i-3 carbonyl oxygen of the backbone could bend the a 
helix creating an angle of 4.3°.Reffff Ballesteros et. al. Ref Ballesteros J.A., Deupi X., 
Olivella M., Haaksma E.E.J, and Pardo L. Serine and threonine residues bend a- helices 
in the %i = g Conformation. Biophysical Journal Volume 79 November 2000-2754-2760 

Presence of proline in a Helices 

It is known that the presence of proline residues in a-helices can perturb the 
structure by introducing a kink between the segments preceding and following the proline 
residues. The distortion of this helical structure is the result of the avoided steric clash 
between the ring of the proline at position (i) and the backbone carbonyl at position (i-4), 
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as well as the elimination of helix backbone H-bonds for the carbonyls at position (i-3) 
and (i-4). Figure 12 shows an example of the structural features of a proline kink. 



N-Terminus 




Bend 
Angle 



Figure 1 1 . The structural features of a Pro kink motif that deviate from a standard a- 
helical conformation are illustrated here. 

The Pro residue (yellow) induces a bend in the a-helix, defined as the angle 
between the axes (yellow) of the helical segments preceding the Pro kink and following 
it. The opening of the helix backbone induced by the Pro residue results in an increased 
solvent accessibility for the carbonyls of amino acids located at positions 1 — 3 and / - 
4 from the Pro. The hydrogen bonding pattern of the a-helix is broken in the Pro kink 
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motif: the carbonyls of amino acids located at positions 1 — 3 and / 
not hydrogen bonded. 

Bend angle and wobble angle definition 



- 4 from the Pro are 



A B 




Figure 12. Illustrates the bend and wobble angle. 



Bend angle: Angle between the two parts when the helix is kinked along its axis. 
Wobble angle: Angle that defines the orientation of the post-proline helix in three 
dimensional space, with respect to the pre -proline helix. 
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FACE SHIFT 

IDEAL OVEft-WOUND UNOER-WOUND 

y y y 




Figure 13. Illustrates the face shift. 

Face shift: Measures the distortion that causes a twist of the helix "face" in such a 
way that amino acids that used to be on the same side (face) of the helix are shifted and 
are on a different side as a results of presence of the proline bend 
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Sequence Dictated Divergence 




CB1 TMH6 DIRLAKTLVLILVVLI ICWGPLLAIMVYDVFG 

CB2 TMH6 DVRLAKTLGLVLAVLLICWFPVLALMAHSLAT 

Figure 14. Illustrates an example of sequence dictated divergences between CB1 and 
CB2 receptor. 

An example of sequence dictated divergences is the study of TMH6 in CB1/CB2 
where the residue 6.49 preceding P6.50 of motif CWXP between CB1 and CB2 is G6.49 
for CB1 and F6.49 for CB2. The conformational memories for CB1 TMH6 showed a 
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wider range of kink angles than CB2. But when the mutants where switch as shown in 
figure 15, CB2 shown the same angles of kink. 

R and R* Model Refinement 

Revised versions of TMHs 2, 5, 6 and 7 that result from CM studies were 
substituted for the existing helix in the GPR55 model. The OPLS_2005 all atom force field 
in Macromodel 9. 1 (Schrodinger, LLC, New York, NY, 2005) was used to energy minimize the 

o 

GPR55 R and R* models. An 8.0 A extended non-bonded cutoff (updated every 10 steps), 

o o 

20.0 A electrostatic cutoff, and 4.0 A hydrogen bond cutoff was used in each stage of the 
calculation and the energy was minimized until an energy gradient of 0.1 kcal/mol is reached. 

Loops affected by the replacement of TMH2 ((IC-1,EC-1), TMH5 (EC-2, IC-3)), 
TMH6(IC-3,EC-3) and TMH7 (EC-3) were modeled using the Modeller program (42,43). 
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CHAPTER IV 
RESULTS 

TMH2 

Conformational memories (CM) calculations on TMH2 of GPR55 were used to 
explore the accessible conformation of TMH2. The average proline kink angle and 
standard deviation for each set of conformational results generated by CM were measured 
using the ProKink program (40). This program, which is part of Simulaid Conversion 
program was used to calculate the face shift, wobble angle, and bend angle of each helix. 
A statistical analysis of the resultant face shift, wobble angle and bend angle for the 105 
helices containing the proline at position 2.58 (78) was performed using the statistical 
analysis package (GraphPad Software for one sample t test) and compared to TMH2 of 
the pVadrenergic receptor which contains a proline at position 2.59. Figure 15 illustrates 
the primary amino acid sequence of TMH2 in GPR55 vs. the pVadrenergic receptor. 



Helix 2 sequence 

2 2 
5 5 
8 9 

GPR55 ATSI YMINLAVFDLLLVLSLPFKMVLSQVQ 
Be t a-2 VTNYFIT SLACADLVMGLAVVPFGAAH I LM 



Figure 15. Primary amino acid sequence of TMH2 in GPR55 and the p2-adrenergic 
receptor. 
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We expected that the presence of proline residues and the different position of this 
residue in the Beta 2 receptor and GPR55 (Beta 2 has proline at 2.59 but GPR55 has 
proline at 2.58) can perturb the helical structure and give a different conformation of 
TMH2. 

We considered the region from V2.54 to P2.58 as flexible helix regions and serine 
2.56 near the Proline 2.58 was set in g- position for Xi torsion angle 
A one sample independent t test was used to analyze the means of the bend angle, wobble 
angle, and face shift for the GPR55 CM results containing P2.58 (78) compared to the 
crystal structure of P2-adrenergic receptor (containing P2.59). This test revealed that, at 
the p=0.05 level, the difference of the population means of GPR55 CM results bend angle 
and face shift were significantly different from TMH2 in the crystal structure of P2- 
adrenergic receptor. No statistically significant difference was found for the wobble angle 
between GPR55 TMH2 CM results and TMH2 in crystal structure of P2-adrenergic 
receptor. Results are shown in Table 3. 



Table 3. Comparison of Pro Kink analysis of CM output vs Prokink analysis of TMH2 in 
the crystal structure of p2-adrenergic receptor 



Proline 
Position 


Bend Angle 
(degrees) 


Face Shift 
(degrees) 


Wobble Angle 
(degrees) 


P2.59 
(P 2 -adrenergic) 


17 


95.4 


-83.3 


*P2.58 (GPR55) 


**39.2 ± 8.5 


**77.9±30.7 


-79.1± 45.4 



*Values are presented as mean + S.D. 

** Statistically significant difference at P 0.05 level in one -sample independent t test. 
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The 105 helices from the CM output were superimposed from A2.38 (58) to 
L2.53 (73) to compare the populations of TMH2s. In the context of a three dimensional 
model of the GPR55 receptor, the criteria used to identify a TMH2 helix for bundle 
incorporation were: 

The helix must conserve the interaction between residue D2.50 (70) and D7.49 
(284) such that a water molecule could connect these two residues (analogous to the key 
interaction observed in Rho). 

The helix must have no van der Waals overlaps with other TM helices in the 
bundle. The extracellular end of TMH2 must be close enough to the extracellular end of 
TMH3 to meet a distance criterion based upon the extracellular loop -1 (EC-1) length. 
The GPR55 EC-1 loop (SPF) contains half of the number of residues as Rho (GYFVFG). 
The attachment of EC-1 loop must be possible without over tightening the TMH 
backbone that would cause van der Waals overlaps. 

The resultant superposition with the CM results superpositioned to the beta-2 
adrenergic receptor is shown in figure 16. 



36 




Figure 16. A side view of the 105 TMH2 helices from the CM runs superimposed on 
TMH2 of the pVadrenergic receptor is illustrated here. These helices were superimposed 
from the residue 2.38 to residue 2.53. 

The strength of the CM method is based in allowing sampling of low free-energy 
conformations of a helix and to assess the ability of a helix-distorting residue (Pro, Gly or 
Ser) to affect the conformation of the helix. The limitation of this method is that the 
calculation is performed on an isolated helix, but the helix needed for the model is one 
that will be part of the TMH bundle and therefore must meet several steric constraints 
imposed by the receptor. Five helices from the CM runs met, in the context of a three 
dimensional model of the GPR55 receptor, the criteria used to identify a helix for bundle 
incorporation. The most important criterion used to select the appropriate helix required 
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that the EC-1 loop attachment be possible without overtightening of the helix ends 
(which would cause van der Waals overlap due to the short three residue EC-1 loop of 
GPR55 (SPF). The helix finally chosen for the bundle incorporation was helix number 40 
from CM with a bend angle of 36.8°, a wobble angle of -81.5°, and a face shift of 62.0°. 
When this helix was incorporated into the TMH bundle, the residues Ser 2.64 (residue 
27), Leu 2.63 and Lys 2.60 were facing into the binding site crevice. 

In helix 2 different GPCRs have a single proline at position 2.58, 2.59 or 2.60. 
Proline can be in different position and act as mimic of GGFT motif present in GPCRs. 
This results showed that shift by one position of the proline modified the angle between 
the pre-proline and post-proline helix. 

TMH5 

CM calculations were used to explore the accessible conformations of GPR55 
TMH5. TMH5 in GPR55 contains the highly conserved P5.50 that is used as an 
alignment guide for most class A GPCR's, but also contains a second proline at position 
P5.41. The average proline kink angle and standard deviation for each set of 
conformational results generated by CM were measured using the ProKink program (40). 
This program, which is part of Simulaid Conversion program was used to calculate the 
face shift, wobble angle, and bend angle of each helix. A statistical analysis of the 
resultant face shift, wobble angle and bend angle for the 105 helices containing the 
proline at position 5.41 (184) was performed using the statistical analysis package 
(GraphPad Software for one sample t test) and compared to TMH5 of the pVadrenergic 
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receptor which contains a proline at position 5.50. See Figure 17 for primary amino acid 
sequence of TMH5 in GPR55 vs. the pVadrenergic receptor. 



Helix 5 sequence 

5 5 
4 5 
1 O 

GPR55 AKVFFP LE VFGF L LPMGI MGFCC SRS I H I LL 
Beta-2 QAYAIASSIVSFYVP LVIMVFVYSRVFQEAK 



Figure 17. Primary amino acid sequence of TMH5 in GPR55 and the pVadrenergic 
receptor. 

As mentioned previously, the presence of proline residues in a-helices can 
perturb the structure by introducing a distortion of this helical structure. This distortion is 
due to the avoided steric clash between the ring of the proline at position (i) and the 
backbone carbonyl at position (i-4), as well as the elimination of helix backbone H-bonds 
for the carbonyls at position (i-3) and (i-4). 

Helix 5 was built based on the backbone of the Beta 2 adrenergic receptor and the 
region around P5.50 was restrained to preserve the initial structure as the Beta 2 
adrenergic receptor. Only the flexible region containing P5.41 was varied because the 
effects of P5.41 on helix geometry are unknown and this can be solved using CM. 
The backbone torsions of residues K5.37 to P5.41 were varied ±50° 

The region surrounding P5.50 in the pVadrenergic receptor is distorted, losing 
helicity in the local region rather than adopting a significant proline kink. This regional 
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distortion was incorporated into GPR55 TMH5 before CM calculations were performed. 
Table 3 lists the proline bend, wobble and face shift angles for the P2-adrenergic receptor 
structure. GPR55 TMH5 has an additional proline residue at 5.41. The region around this 
proline was varied during CM calculations (see Methods). Table 4 lists the values for the 
bend, wobble and face shift angles for GPR55 TMH5 in the region of P5.41 . It is clear 
here that P5.41 induces a different helix distortion compared to P5.50 in the P2- 
adrenergic receptor. A comparison of the distortions is provided in Table 4 by using a one 
sample independent t test to analyze the means of the bend angle, wobble angle, and face 
shift for the GPR55 CM results containing P5.41 (184) compared to the crystal structure 
of P2-adrenergic receptor (containing P 5.50). This test revealed that, at the p=0.05 level, 
the difference of the population means of GPR55 CM results bend angle, wobble angle 
and face shift were significantly different from the TMH5 P5.50 region in the crystal 
structure of P2-adrenergic receptor. 



Table 4. Comparison of ProKink analysis of CM GPR55 TMH5 output vs Prokink 
analysis of TMH5 in the crystal structure of p2-adrenergic receptor. 

Proline Position (degrees)^ 6 ^ aCC (degreed Wobble Angle 




P5.50 (211) 
(P2-adrenergic) 


6.0 


79.0 


-137.2 


*P5.41 (184) (GPR55) 


**23.2 ±7.3 


**34.9± 16.5 


-67.1+ 138.9 



* Values are presented as mean ± S.D. 

** Statistically significant difference at P 0.05 level in one-sample independent t test. 



40 



Figure 18. Illustrated here are side and extracellular views of the 105 TMH5 helices from 
the CM runs superimposed on TMH5 of the pVadrenergic receptor (from the residue 5.53 
to residue 5.64). 

Rho and p2-adrenergic contain the conseved proline residue at position 5.50 
compared to GPR55 which contains two proline residues. The conformational memories 
results showed a distortion of the helical structure to avoid steric clash between ring of 
proline and backbone carbonyl, contributing to the observed flexibility. 

TMH6 

TMH6 is an important helix in the R to R* transition, for this reason I studied the 
conformations accessible to TMH6 in GPR55 using CM. TMH6 in GPR55 does not have 
the class A GPCR highly conserved CWXP motif , but it contains the conservative 
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substitution in this motif: SFLP. See Figure 19 for the primary amino acid sequence of 
GPR55 TMH6 vs. that of the P2-adrenergic receptor. 



Helix 6 sequence 

6 6 
4 5 

7 O 

GPR55 QKACI YSI AASLAVFVVSFLPVHLGFFLQFLV 
Beta-2 EHKALKTLG 1 1 MGTFTLCWLPFFI VN I VHVIQ 



Figure 19. Primary amino acid sequence of TMH6 in GPR55 and the P2-adrenergic 
receptor. 

The average proline kink angle and standard deviation for each set of 
conformational results generated by CM were measured using the Pro Kink program (40). 
This program, which is part of Simulaid Conversion program was used to calculate the 
face shift, wobble angle, and bend angle of each helix. A one sample independent t test 
was used to analyze the means of the bend angle, wobble angle, and face shift of the 105 
helices containing proline at position 6.50 (GPR55 P6.50 (241)) and compared to the 
crystal structure of P2-adrenergic receptor ( P2-adrenergic P6.50 (288)). This test revealed 
that, at the p=0.05 level, the difference of the population means of the GPR55 CM results 
for wobble angle and face shift were significantly different from TMH6 in the crystal 
structure of P2-adrenergic receptor, however, the bend angle was not significantly 
different (see Table 5). 
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Table 5. Comparison of Pro Kink analysis of Conformational Memories output vs Prokink 
analysis of TMH6 in the crystal structure of p2-adrenergic receptor. 



Proline Position 


(degrees) 


(degrees) 


(degrees) 


P 6.50 (288) 










33.2 


84.9 


-73.0 


(P2-adrenergic) 








*P6.50 (241) 
(GPR55) 


32.6 ± 11.4 


55.7 + 38.7 


-90.1 +74.2 



Values are presented as mean + S.D. 

** Statistically significant difference at P 0.05 level in one-sample independent t test. 




Figure 20. Illustrated here are side and extracellular views of the 105 TMH6 helices from 
CM runs for the GPR55 inactive state superimposed on TMH6 of the P2-adrenergic 
receptor (from the residue 6.32 (223) to residue 6.45(236)). 

In the P2-adrenergic receptor, aromatic residues at 6.44 and 6.52 flank the toggle 

switch residue, W6.48. These residues, along with C6.47 have been proposed to be part 
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of a larger rotamer toggle switch within the binding site crevice of the P 2- adrenergic 
receptor (39). GPR55 also has aromatic residues flanking F6.48, F6.44 and H6.52, while 
the residue at 6.47 is a serine. To test the relation of the Serine 6.47 with the modulation 
of a possible GPR55 rotamer toggle switch (F6.44/F6.48/H6.52), the population of the %l 
rotamer for each starting conformation based on the position of S6.47 was correlated with 
the position of the residues F6.48, F6.44 and H6.52 (see Table 6). 

Table 6. The correlation of %1 of S6.47, F6.44, F6.48 and H6.52. 



Ending 


F6.48 


F6.44 


H6.52 


rotamer 


g- 


g+ 


t 


g- 


g+ 


t 


g- 


g+ 


t 


S6.47 


g+ 


0% 


26% 


74% 


0% 


15% 


85% 


0% 


25% 


75% 


t 


0% 


26% 


74% 


0% 


10% 


90% 


0% 


22% 


78% 


Javitch 


W6.48 


F6.44 


F6.52 


results 


g- 


g+ 


t 


g- 


g+ 


t 


g- 


g+ 


t 


C6.47 


g+ 




39% 


61% 










23% 


77% 


t 




70% 


30% 










44% 


56% 



The results were compared to those reported by Javitch and Ballesteros in 2002 
for P2 -adrenergic receptor (P2-AR) (39). These results indicate that there is no correlation 
between the xl of S6.47 and the xl of the residues F6.44 (235), F6.48 (239) and H6.52 
(243). 

The next step was to test if there is a correlation between the xl of F6.48 and the 
Xl's of F6. 44(235) and H6. 52(243). The results of this study are summarized in Table 7. 
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Table 7. The correlation of y\ of F6.48 (239) with F6.44 (235) and H6.52 (243) 



y\ of serine 6.47 (238) was held in trans position 




F6.44 


H6.52 




gplus 


trans 


gplus 


trans 


Bend angle 


Wobble angle 


F6.48 


gplus 


30.% 


70.% 


81.% 


19.% 


*40.5 ± 16.3 


*-56.7 ± 60.5 




trans 


4.% 


96.% 


3.% 


97.% 


*29.9 + 7.6 


*-101.7±74.7 



*Values are presented as mean ± SD. 



The results indicate that there is a correlation between the xl of F6.48 (239) and 
%\ of H6.52. When the %1 of F6.48 (239) is in g-plus (inactive state conformation), the 
Xl of H6.52 prefers the g plus position. When the xl of F6.48 (239) is in trans position 
(activated state), the xl of H6.52 prefers the trans position. Results for F6.44, suggest 
that trans position is favored for F6.44 whether the xl of F 6.48(239) is g-plus or trans. 
There is, however, a lessening of the number of helices with xl of F6.44 in trans when 
the xl of F6.48 is g-plus (96% 70%), but a trans xl for F6.44 is clearly still favored. 

Results in Table 7 also show that although there is a trend in the bend angle of 
TMH6 towards less bent helices when F6.48 undergoes a xl g+— ► trans change, this is 
not statistically significant. This result is not consistent with experimental evidence that 
TMH6 straightens during activation. 

TMH7 

Figure 21 illustrates the sequence alignment for GPR55 TMH7 vs. that of the pV 

adrenergic receptor. TMH7 in Class A GPCRs usually has the conserved sequence motif, 

NPXXY, which not only influences the conformation of TMH7, but also places Y7.53 in 

the correct position to interact with F7.60 on the intracellular extension of TMH7, Hx8. 
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The preferred conformation of TMH7 in GPR55 was studied because it lacks the 
NPXXY motif in TMH7. Instead, this region of the sequence is DVFCY in GPR55. 
C7.46 to V7.50 was considered the flexible region in CM calculations. The output of the 
CM was superimposed on the P2 -adrenergic receptor (pVAR) template from residue 
F7.51 (286) to V7.56 (291) (intracellular region) (Figure 21). 



Helix 7 sequence 

7 
5 


GPR55 FLQLSMCFSNVNCCLDVFCYYFVI 
Beta-2 VY I LLNWI GYVNSG FNPLI YCRD P 



Figure 21. Primary amino acid sequence of TMH7 in GPR55 and the pVadrenergic 
receptor. 
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Figure 22. Side view of the TMH7 CM output helices superimposed to the (pV 
adrenergic receptor. 



Table 8. Comparison of Pro Kink analysis of CM output vs Prokink analysis of TMH7 in 
the crystal structure of p2-adrenergic receptor 



Proline Position 


(degrees) 


(degrees) 


(degrees) 


P7.50 

(P2-adrenergic) 


20.1 


9.1 


161.4 


*V 7.50 (GPR55) 


**10.1 +4.7 


16.9+144.2 


**16.8+ 11.6 



* Values are presented as mean + S.D. 

** Statistically significant difference at P 0.05 level in one-sample independent t test. 

[3-2 adrenegic receptor contains a proline in position P7.50 which forms part of 
the motif NPXXY but GPR55 lacks the proline residue in position P7.50. The output of 
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the conformational memories showed a decreased bend compared to the [3-2 adrenegic 
receptor. 

Helix 8 GPR55 

The intracellular extension of TMH7, Hx8 ((E 7.59 (294) - A 7.86 (319)) was 
modeled in the same conformation seen is the Rho structure and the receptor was 
truncated after A7.86. 

Model of Inactive state (R) and active state form of GPR55 

A model of the R and R* form of GPR55 was created using the crystal structure 
of the pVadrenergic receptor as a starting point. The sequence of the human GPR55 
receptor (see Figure 24) was aligned with the sequence of the pVadrenergic receptor 
using highly conserved residues across Class A GPCRs as alignment guides. Results 
from CM were used to substitute TMH2, TMH5, TMH6 and TMH7 for existing helices 
in GPR55 model. Loops affected by the replacement of TMH2 (IC-1, EC-1), TMH5 (EC- 
2, IC-3), TMH6 (IC-3, EC-3), TMH7 (EC-3) were modeled using Modeller program. 
OPLS_2005 all atom force field in Macromodel 9.1 was used to energy minimize the 
GPR55 R and R* models. 
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Intracellular 



Figure 23. Schematic representation of human GPR55 receptor. The most highly 
conserved residue in each helix is highlighted in blue. 

The biophysical studies reporting the changes that occur during the R to R* 
transition from Rho and the P2 -adrenergic receptor are the bases for the creation of the 
GPR55 model active state (R*) (19). These studies indicated the rotation of TMH3 and 
TMH6 and a conformational change in TMH6 (straightening in the CWXP hinge region 
of TMH6) which occurs upon activation. An R* GPR55 model was created including a 
TMH6 conformer derived from our conformational memories study of GPR55 TMH6. 
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CHAPTER V 
DISCUSSION 



In the present study, I used the Conformational Memories (CM) technique to 
explore the sequence dictated conformation of key helices: TMH2, TMH5, TMH6 and 
TMH7 in GPR55 and compared these results to the p 2 -adrenergic receptor structure (44). 

TMH2 

Many GPCR's contain Pro at position 2.59 or 2.58 so the bend is present at the 
residues 2.55 to 2.59 or 2.58 in the helix turn formed by the Pro-kink at the similar 
positions as the bend in rhodopsin, which contains the motif Gly-Gly-X-Thr-Thr instead 
of proline. 

Different GPCR's have a single Pro at position 2.58, 2.59 or 2.60 and each one of 
these positions would be expected to orient the bend differently. The fact that the Pro can 
be in different positions and act as mimic of the motif Gly-Gly-X-Thr-Thr suggest that 
this may be a source of considerable structural diversity in the Class A receptors. 

Our CM analysis of TMH2 in GPR55 compared to the TMH2 of the pY 
adrenergic receptor showed that the shift by one position of the amino acid proline 
(GPR55 containing P2.58 (78) compared to the crystal structure of pVadrenergic receptor 
containing P2.59) produces a significant difference in the bend angle and face shift of 
TMH2. 
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Meaning that the different position of proline modified the angle between the two 
parts (pre -proline and post-proline) of the TMH2 in GPR55 compared to the pV 
adrenergic receptor, causing a shift of the amino acids as a result of the presence of the 
proline bend. The presence of the serine residue in position 2.56(76), satisfied its 
hydrogen-bonding capacity by an intrahelical hydrogen bond interaction between the O-y 
atom and the i-4 carbonyl oxygen of the residue Leu 2.52(72) stabilizing the bend angle 
created by the presence of the proline. 

TMH5 

In Rho and the pVadrenergic receptor there is a single proline residue in TMH5 at 
position 5.50, which causes a deviation from normal alpha-helicity in the region, 5.46- 
5.50. TMH5 in GPR55 has two proline residues (P5.40 and P5.50). The CM results 
showed a difference of the population means of GPR55 in the bend angle, wobble angle 
and face shift (for P5.40) compared to the TMH5 in the crystal structure of pVadrenergic 
receptor(for P5.50). The distortion of the helical structure results from the avoidance of a 
steric clash between the ring of the proline at position (z) and the backbone carbonyl at 
position (z - 4), as well as the elimination of helix backbone hydrogen bonds for the 
carbonyls at positions (i - 3) and (z - 4). These two factors (departure from the ideal 
helical pattern and the reduction in hydrogen-bond stabilization) contribute to the 
observed flexibility of a proline-containing -helix, which is illustrated in Figure 24. 
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Figure 24. Presence of proline causes a distortion of the helical structure. This figure 
shows how the presence of proline at position 5.50 causes the helical portion to be 
disrupted between the leucine 5.48 (#13 in figure) and phenylalanine 5.47 (#12 in figure). 

TMH6 

TMH6 in GPR55 does not have the Class A GPCR highly conserved CWXP 
motif. However, it has the conservative substitutions of this motif: SFLP. Since TMH6 
has been implicated in conformational changes during receptor activation, it was 
important to study the allowed conformations of TMH6. It has been reported during the 
activation process of the GPCRs, that agonist binding influences the rotameric state of 
aromatic residues in TM6 altering the configuration of TMH6 prokink which in turn 
results in movement of the cytoplasmic end of TMH6 away from TMH3. For the pV 
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adrenergic receptor, the residue C6.47 has been reported to modulate the configuration of 

aromatic cluster surrounding the proline kink and the proline kink itself (45). The cluster 

of highly conserved aromatic residues surrounding the proline kink includes Phe 6.44, 

Trp 6.48 and Phe 6.52 in the P2 -adrenergic receptor, all of which face into the binding 

site crevice. In the case of GPR55, the residues forming the aromatic cluster are F6.48, 

F 6.44 and H 6.52. To test the ability of Serine 6.47 to modulate the aromatic cluster (as 

suggested by the results of Shi and co-workers for the p 2 -adrenergic receptor (Shi et al., 

2002)), the population of the xl rotamers for each starting conformation based upon the 

rotameric state of S 6.47 was correlated with the position of the residues F 6.48, F 6.44 

and H 6.52 (see Table 5). These results did not show a correlation between the position of 

the rotameric state of S6.47 and the %1 position of the residues F6.44, F 6.48 and H 6.52. 

However, there was a correlation between the rotameric state of residue H6.52 based on 

the position of F 6.48 (See results in Table 6). When the xl of F6.48 (239) is in g-plus 

(inactive state conformation), the xl of H6.52 prefers the g plus position. When the xl of 

F6.48 (239) is in trans position (activated state), the xl of H6.52 prefers the trans 

position. In contrast, results for F6.44, suggest that trans position is favored for F6.44 

whether the xl of F 6.48(239) is g-plus or trans. The resulting simulation structure can be 

clustered into two major groups based in the rotamer conformation of 6.48 and 6.52: (1) 

F 6.48 g plus/ H 6.52 g plus and (2) F 6.48 trans/H 6.52 trans. 

Although the results in Table 6 suggest a trend in the F 6.48 xl g+ — > trans 

change towards less bent helices, this trend is not statistically significant. This result, 

therefore, is not consistent with biophysical studies that show that during activation 
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P6.50, in the region CWXP region of TMH6 can act as a flexible hinge allowing TMH6 
to straighten upon activation. This movement would result in its intracellular end moving 
away from TMH3 and upwards towards the lipid bilayer. 




Figure 25. Ionic lock in GPR55 inactive state 



There is evidence that activation is accompanied by rigid domain motion of 
TMH3 and TMH6. In the intracellular ends of TMH3 and TMH6 there is a constrained 
by D3.49 (118); R3.50 (1 19) and Q6.30 (221) forming a salt bridge that limits the 
mobility of the cytoplamic ends and acts like ionic lock. During activation, P6.50 acts as 
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a flexible hinge, permitting TMH6 to straighten upon activation and move its intracellular 
end away from TMH3 and up toward the lipid bilayer 

The residue F6.48 has been reported to change conformation upon activation: 
R y\ gn — > trans in R* 

Many GPCRs have aromatic residues that flank W6.48. acting as a toggle switch 
with possible correlated movements of flanking aromatic residues such as F6.44 and 
H6.52 in GPR55 as shown in figure 26. 




Figure 26. Toggle switch in GPR55. 
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The correlation of xi of F6.48 (239) with F6.44(235) and H6.52 (243) are reported 
in table 9. 

Table 9. Correlation of Xi of F6.48 (239) with F6.44(235) and H6.52 (243). 



%1 of serine 6.47 (238) was held in trans position 




F6.44 


H6.52 




gplus 


trans 


gplus 


trans 


Bend angle 


Wobble angle 


F 

6.48 


g 

plus 


30% 


70% 


81% 


19% 


*40.53± 16.26 


*-56.73± 60.50 




trans 


4% 


96% 


3% 


97% 


*29.91 ±7.64 


*-101.65±74.67 



Values are presented as mean + SD 



The results shown: 

• F6.44 prefers trans: Population does not correlate with F6.48 %1 

• H6.52 %l correlates with F6.48 %l 

• When F6.48 xl is g+, H6.52 xl prefers g+ 

• When F6.48 xl is trans, H6.52 xl prefers trans 

The main point used to assemble de model of active and inactive GPR55 receptor 
were: 

• No Van der Waals overlaps with other TMHs in the bundle 

• The extracellular end of TMH2 must be close enough to the extracellular end of 
TMH3 

• The GPR55 EC-1 loop (SPF) contains half of the number of residues as Rho 
(GYFVFG). 
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• The attachment of EC-1 loop must be possible without over tightening the TMH 
backbone that would cause Van der Waals overlaps. 

• Based on biophysical literature on the R to R* transition in Class A GPCRs, 

• Activation is accompanied by a straightening of TMH 6 

• Activated model of GPR55 

• conformer chosen for R* bundle in which the proline kink angle is moderated in 
TMH 6 

• F6.48 xl is in trans and H6.52 %1 is in trans 

Following these points we obtained the models shown in figure 27. 




Figure 27 Inactive receptor GPR55 top view ( left figure) and active receptor GPR55 
(right figure). 

TMH7 

The preferred conformation of TMH7 in GPR55 was studied because it lacks the 
NPXXY motif in TMH7. Instead, this region of the sequence is DVFCY in GPR55. The 
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output of the CM showed that without the presence of proline, the helix is less bent 
compared to the beta adrenergic receptor TMH7. 
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CHAPTER VI 
CONCLUSIONS 



In conclusion, the results of the CM analysis of GPR55 are consistent with the 
proposal of Ballesteros et al.(46) , that even though the overall structure of Rho and other 
class A GPCRs may be very similar, there are localized regions where the structure of 
these receptors diverges. A significant range of conformational diversity could be 
generated within the binding site crevice of different GPCRs by the presence of Pro-kinks 
and Cys/Ser/Thr residues. The presence of different amino acids or alternate 
microdomains may support very similar deviations from regular a-helical structure, 
resulting in a similar tertiary structure. This could be a mechanism by which the different 
receptors could diverge sufficiently to develop the selectivity necessary to interact with 
diverse ligands, while still maintaining a similar overall fold in the core function. The 
core function is the signaling activation through a conformational change in the TMHs 
that alters the conformation of the cytoplasmic surface and subsequent interaction with 
G-proteins. This core function is shared by the entire Class A GPCR family of receptors, 
despite their selectivity for a diverse group of ligands. In this thesis, the consequence of 
sequence divergences in GPR55 was explored. The results reported here should help to 
define the mechanism of drug-receptor interaction relevant to cannabinoid physiological 
and pathophysiological functions including drug abuse. 
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